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ABSTRACT 

A prevalent but untested paradigm is often used to describe the prominence-cavity system: the cavity is under-dense 
•because it is evacuated by supplying mass to the condensed prominence. The thermal non-equilibrium (TNE) model 
of prominence formation offers a theoretical framework to predict the thermodynamic evolution of the prominence 
and the surrounding corona. We examine the evidence for a prominence-cavity connection by comparing the TNE 
model with diagnostics of dynamic extreme ultraviolet (EUV) emission surrounding the prominence, specifically 
prominence horns. Homs are correlated extensions of prominence plasma and coronal plasma which appear to 
connect the prominence and cavity. The TNE model predicts that large-scale brightenings will occur in the Solar 
Dynamics Observatory Atmospheric Imaging Assembly 171 A bandpass near the prominence that are associated 
with the cooling phase of condensation formation. In our simulations, variations in the magnitude of footpoint 
heating lead to variations in the duration, spatial scale, and temporal offset between emission enhancements in 
the other EUV bandpasses. While these predictions match well a subset of the horn observations, the range 
of variations in the observed structures is not captured by the model. We discuss the implications of our one- 
dimensional loop simulations for the three-dimensional time-averaged equilibrium in the prominence and the 
cavity. Evidence suggests that horns are likely caused by condensing prominence plasma, but the larger question of 
whether this process produces a density-depleted cavity requires a more tightly constrained model of heating and 
better knowledge of the associated magnetic structure. 
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1. INTRODUCTION 

Prominences are high-density, low-temperature structures 

that are suspended in the corona. They are often surrounded by a 

low-density coronal structure known as a cavity (Saito & Hyder 

1968). Two fundamental unsolved problems in prominence 
physics are: how is the prominence supported and how is 
mass supplied to the prominence? The commonly accepted 
answer to the first question is that the magnetic field supports 
the prominence mass and must be dipped to collect stable, 
condensed plasma. Two models were proposed for the geometry 
of dipped fields: sheared arcades (Kippenhahn & Schliiter 1957; 
Antiochos et al. 1994) or flux ropes (Kuperus & Raadu 1974; 
Low & Hundhausen 1995). The flux-rope model may also 
explain the morphology of the cavity — helical field lines on the 
interior of the flux rope form an elliptical structure when viewed 
down the axis of the flux rope (Gibson & Fan 2006). The three- 
dimensional (3D) magnetic geometry of the prominence-cavity 
system has been addressed through MHD (Fan & Gibson 
2003 ; DeVore & Antiochos 2000) and magnetofrictional models 
(Mackay & van Ballegooijen 2001), neither of which has 
incorporated an accurate coronal energy treatment thus far. 

The question of how mass is supplied to the prominence has 
also been theoretically addressed. The thermal instability of 
coronal energy balance suggests that coronal equilibrium is sus- 
ceptible to large-wavelength, isobaric condensing perturbations 
(Field 1965). At present, the thermal non-equilibrium (TNE) 
model represents the most detailed theoretical explanation of 
the formation of a prominence through a form of thermal in- 
stability. This model predicts that plasma condensations form 
through time-dependent evaporation and cooling driven by heat- 


ing localized at loop foolpoints (Mok et al. 1990; Antiochos & 
Klimchuk 1991; Kuin & Martens 1982). The geometry of each 
flux tube dictates the properties of each condensation, including 
whether condensations form and remain stationary or undergo a 
cycle of formation and destruction by falling onto the chromo- 
sphere (Karpen et al. 2001, 2005), while the heating timescale 
determines whether condensations form at all (Karpen & 
Antiochos 2008). TNE in prominences, coronal loops, and coro- 
nal rain has been studied in detail using one-dimensional (ID) 
hydrodynamic and multidimensional MHD simulations (Karpen 
et al. 2003, 2006; Muller et al. 2003; Mok et al. 2008; Luna et al. 
2012). The dynamic evolution of density and temperature de- 
termines the emissivity of the plasma along the loop, which 
in turn is directly comparable to observations through forward 
modeling. 

Many observational studies have constrained the plasma 
properties within the cavity as well as the prominence-corona 
transition region (PCTR). The large-scale temperature and 
density structure of the cavity has been documented by Schmit 
& Gibson (2011) and Kucera et al. (2012). Emission inside the 
cavity from ion species characteristic of hot plasma (relative 
to quiet Sun) have been observed by Habbal et al. (2010) 
and Reeves et al. (2012). The dynamics and fine structure of 
the PCTR have been discussed by Heinzel & Anzer (2001) 
and Heinzel et al. (2008), in which their analysis focuses 
on hydrogen and helium lines and continua forming between 
the 10 4 K condensation and the 10 6 K corona. A complete 
review of spectral diagnostics on prominences and their coronal 
environments is found in Labrosse et al. (2010). 

In Schmit & Gibson (2013), hereafter referred to as Paper 1, 
we characterized the properties of time-dependent extreme 
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Figure 1, Overlay of 171 A horns over a 193 A image of the cavity observed 
2011 February 28. A detailed discussion of the extraction and analysis of horn 
structures is presented in Paper 1 . 

(A color version of this figure is available in the online journal.) 

ultraviolet (EUV) emission in a prominence and the asso- 
ciated cavity observed by the Solar Dynamics Observatory 
Atmospheric Imaging Assembly (SDO/ AIA) (Lemen et al. 
2012). In this research, we compare these observations with 
emission structures predicted by the TNE model to derive insight 
into the poorly understood relationship between the prominence 
and cavity. 

In Section 2, we briefly summarize the observation of promi- 
nence horns. In Section 3, we describe the hydrodynamic 
simulation, the emission structures that are produced during 
the TNE process, and how that emission compares with horn 
observations. In Section 4, we consider the roles of horns and 
the TNE process in producing the low-density cavity. Our con- 
clusions and the next steps necessary to understand the dynamic 
prominence-corona system are summarized in Section 5. 

2. HORN OBSERVATIONS 

Paper 1 details the diagnostic methods used to characterize 
prominence horns from multi-bandpass SDO/AIA data sets, 

so we summarize those results briefly here. The dominant 

components of emission variability in the prominence-cavity 
region are prominence horns (Paper 1; Section 2). Figure 1 
shows an overlay of the outlines of several horns extracted from 
a cavity on 201 1 February 28. Homs are most clearly visible 
in the 171 A bandpass as strong brightenings that extend from 
the prominence into the cavity, achieving heights similar to the 
50 Mm coronal hydrostatic scale height. The 171 A brightenings 
were correlated with changes in the perimeter of the prominence, 
such that 304 A emission would appear to extend up the base 
of the horn, and were accompanied by weak 193 A bandpass 
brightening. Temporally, we found that that the peak brightness 
of the 171 A emission preceded the peak extension of the 304 A 
prominence by 30 ± 18 minutes, whereas a —0.6 d= 23 minute 
time lag existed between 193 A and 171 A emission. 

In addition to temporal and spectral information, prominence 
horns were found to have a particular geometry. A majority 
of prominence horns exhibited a concave-up geometry, which 
appeared constant over the duration of the feature. Because 
these features are likely to be low-/} field-aligned structures, 
they can be modeled with ID hydrodynamic simulations of 
plasma within a coronal flux tube. To present the data in a 
manner that is most easily comparable to our loop models, we 
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created time-distance diagrams similar to those in Karpen et al. 
(2005) and Berger et al. (2012) where the temporal evolution 
of emission along a curved slit is plotted against distance 
along the slit. Figure 2 shows two examples of the 45 horns 
that were analyzed in Paper 1, illustrating their fundamental 
characteristics and the variations between individual horns. 
Figures 2(a) and (c) show the correlation between extensions 
of the prominence in 304 A and larger brightenings in the 
171 A bandpass. The durations of the 171 A brightenings for 
these two horns are nearly identical despite their different 
lengths. Both brightenings occur, simultaneously across the 
entire horn, so there is not a clear prominence-to-cavity or 
cavity-to-prominence progression. However, these horns differ 
significantly in the relationship between the 171 A and 304 A 
emissions. For Horn 10, the 304 A emission reaches maximum 
extension 0.2 hr before the peak emission in 171 A, while for 
Horn 3, the 304 A emission goes through two distinct periods of 
motion: 0 ^ t ^ 0.6 hr and 1 hr ^ t ^ 2 hr. The peak extension 
(which lasts from t = 1.2 hr to £ = 2 hr) occurs 0.6 hr after the 
peak emission in 171 A. This variation in timing is reflected in 
the temporal error bars presented in Paper 1 (Section 3.2). 

Figures 2(b) and (d) are time-distance diagrams for 193 A 
emission, also displaying the 304 A emission in red hatching. In 
both horns, the peak 193 A emission is significantly dimmer 
than the 171 A emission, and the position of peak 193 A 
brightness sits higher along, the slit than the 171 A feature. 
In Horn 10, the peak 193 A emission occurs approximately 
0.1 hr before the 304 A peak extension (171 A and 193 A peaks 
are simultaneous). In Horn 3, the peak 193 A emission occurs 
0.2 hr after the peak 304 A extension and 1 .8 hr after the 171 A 
peak. Thus, the two horns displayed in Figure 2 exhibit different 
correlations between the cooler 171 A emission and the hotter 
193. A emission. 

In summary, we have found that the properties of the 171 A 
emission in horns are roughly uniform throughout the data sets: 
peak intensity and duration do not vary greatly among horns 
despite differences in length, and the emission does not exhibit 
a preferred propagation direction. However, we find significant 
variation in the correlations among the 171 A, 304 A, and 193 A 
emissions. Intense brightenings in the 171 A are unique to 
that bandpass and the properties of that emission are largely 
uniform. The other bandpasses exhibit a range of temporal 
associations with the 171 A feature but a more consistent spatial 
association. As we compare these observations to models of 
the hydrodynamic processes, the degree of temporal correlation 
among these varying emissions will be a primary constraint. 

3. THERMAL NON-EQUILIBRIUM AND 
PROMINENCE HORNS 

The corona is dominated by the magnetic field and the 
timescales for its non-eruptive evolution ranges from days to 
weeks. The thermodynamic evolution time is very rapid in com- 
parison, as thermal conductivity is very efficient at transport- 
ing energy along the magnetic field and the coronal radiative 
loss time is also relatively fast. Assuming a constant density 
of 10^ cm -3 and a radiative loss rate of 10“ 21 erg s" 1 cm 2 3 * , it 
would take approximately 2000 s to cool 10 6 K plasma to 10 5 K. 
Thus, coronal loop models that solve for the evolution of plasma 
contained within a fixed magnetic geometry are appropriate for 
studying the observed prominence and cavity plasmas. 

As in prior investigations of prominence mass formation 
and evolution (e.g., Karpen et al. 2006; Luna et al. 2012), 
we have used the ARGOS code (Antiochos et al. 1999) and 
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Figure 2. Two extracted horn data sets presented in a time-distance plot The x axis is position along the slit (see Paper 1, Section 3.1), and the y axis is time. Panels 
(a) and (c) show 171 A emission in grayscale, (b) and (d) show 193 A emission in grayscale. The red-hatched region denotes where strong 304 A emission is observed. 
(A color version of this figure is available in the online journal.) 


associated post-processing routines to study the dynamics of a 
catastrophically cooling loop and to synthesize its appearance in 
EUV emission lines. In contrast to earlier TNE studies, however, 
the present work is focused on comparing the predicted emission 
with multi wavelength observations of prominences, prominence 
horns, and the surrounding cavities, as presented in Paper 1. 

3.1. Model Initial Conditions 

For purposes of this comparison, we have extracted a deeply 
dipped field line from a stable equilibrium of a partially emerged 
3D flux rope (Gibson & Fan 2006), which is shown in Figure 3. 
Condensed plasma can collect and maintain a quasi-stable 
position within the deep magnetic dip. Most of the dipped 
field line geometries from this flux rope are asymmetric, with a 
high-arched side and a low-flat side. Although those magnetic 
geometries are more prevalent in the 3D flux-rope model, we 
have used a nearly symmetric loop for simplify the analysis of 
the loop emission. 

The ARGOS code solves the following hydrodynamic 
equations: 

T t P + Tt fV = 0 (1) 

dt 3 s 

^(pV) + oV 2 + p ) = -pgiUt) (2) 

dt os 

3 3 / f- - n 3 \ 3 

— U + —IUV- 10~ 6 T 5/2 —T = -P—V + E(t, s) 

dt 3s \ 3s / 3s 

-« 2 A(D, (3) 

where s is distance along the field line, t is time, p is mass density, 
V is the velocity, P is the gas pressure, g|| is the component of 
gravity parallel to the slope of the field line, U is the internal 


energy, E is the heating function, and A is the radiative loss 
function (see Antiochos et al. 1999, for details). The loop code 
uses the radiative loss function of Klimchuk & Cargill (2001). 
This loss function was also used in Karpen et al. (2005), and 
those authors found that form of radiative loss function can 
impact the rate of condensation. While ARGOS is capable of 
including the effects of variations in the cross-sectional area 
of the flux tube, for simplicity in this work we have assumed 
the flux tube has a uniform cross section. Previous studies have 
shown that variations in cross section may affect the rate of 
condensation (Lionello et al. 2013). 

Our dipped loop is initialized with a uniform heating rate of 
6 x 10“° erg cm” 3 * s -1 , and the calculation is continued for 
10 5 s until hydrostatic equilibrium is reached. Over the next 
1000 s, the heating is linearly changed to a combination of 
strong footpoint heating E\ and weak uniform heating Eq in the 
form 

E(s) = E 0 + Ei exp [s t /kl (4) 

where Eq is 6 x 10 -7 erg cm -3 s -1 , s t is the position relative 
to the nearest footpoint, and A is 25 Mm, This choice of A is 
larger than that used in previous experiments due the longer 
loop used for the model. The ratio of length-to-heating scale 
is 21, which comparable to Kaipen et al. (2001; ratio of 22) 
and Mtiller et al. (2004; ratios between 8 and 50). Very short 
scale lengths can lead to formation of multiple condensations 
(see Muller et al. 2004, Figure 13), which would complicate the 
emission analysis. 

By initiating with a high uniform heating, we ensure that the 
initial relaxed loop achieves coronal temperatures. The present 
study is concerned only with the plasma characteristics during 
and after condensation, not with the buildup period. Therefore, 
when the footpoint heating is applied, the background heating 
rate Eq is lowered by an order of magnitude compared with 
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Figure 3. Geometry of a ID symmetric, dipped field line used for the TNE loop simulations, (a) Representative field lines (black) and photospheric magnetic field 
(coior surface) from the Gibson & Fan (2006) flux rope; the extracted field line is shown in red. (b) Height as a function of position along the extracted field line. 

(A color version of this figure is available in the online journal.) 


the initial equilibrium value in order to obtain a condensation 
as rapidly as possible. The lower value of E 0 would produce 
a loop center density and temperature of 2.3 x 10 7 cm -3 and 
0.87 MK, respectively, with, an integrated uniform heat input of 
3 x 10 4 erg cm“ 2 s“* (well below typical coronal values). 

In our analysis, the loop geometry and the basic form of the 
heating function remain unchanged so we can focus primarily 
on the effects of varying E\. We performed two model runs: 
Run A uses E\ = 1.6 x 10~ 3 erg cm -3 s -1 while Run B 
uses Ei — 2 x 10 -4 erg cm -3 s -1 . Runs A and B thus 
have energy inputs of 8 x 10 6 and 1 x 10 6 erg cm -2 s" 1 , 
respectively, comparable to standard coronal heating estimates 
(Withbroe & Noyes 1977). In both cases the footpoint heating 
term dominates the integrated heat input. The variation in 
footpoint heating magnitude between the models results in 
different thermodynamic loop evolution, which affects the 
predicted emission structures . that are compared with hom 
observations. We have examined a range of E\ that produces 
post-condensation coronal plasmas spanning the temperature 
range between the peak temperature responses of - the AIA 
171 A and 193 A bandpasses. This allows us to map the relative 
response in these two bandpasses for the full range of coronal 
loop temperatures that could be present in the quiescent cavity 
(0.8-2.0 MK). Much cooler plasmas are present during the 
dynamic cooling phase. 

3.2. Simulation Results 

■The model coronal loop evolves through three phases, illus- 
trated in Figure 4. The footpoint heating is turned on at t = 0 
(Figure 4, black curves), causing a rapid temperature increase 
that is proportional to E\k. Within 1-2 hr. Run A reaches a peak 
temperature of 3.4 MK (orange curves), while Run B reaches 
1.7 MK. The second phase (purple curves) is characterized by 
increasing density in the coronal segment of the loop evapo- 
rated from the chromosphere, accompanied by increased radia- 
tive losses and large-scale (loop-wide) cooling. This cooling 
phase ends after ~4 — 10 hr with a rapid influx of plasma into 
the loop midpoint, forming a condensation. During the third 
phase (green curves), the loop reaches a new equilibrium with 
hot coronal plasma on both sides of the condensation; for the 
dipped geometry presented in Figure 3 and.the symmetric heat- 
ing profile of Equation (4), the condensation is gravitationally 
stable but continuously growing. At t = 8 hr, the loop resem- 
bles two separate coronal loops. In addition to the usual two 



Figure 4. Density and temperature along the field line at four representative 
times in Run A: initialization (black line), temperature maximum (orange), 
cooling (purple), and post-condensation equilibrium (green). 

(A color version of this figure is available in the online journal.) 


chromosphere-corona transition regions located at the magnetic 
footpoints, each loop has PCTRs located on either side of the 
condensation. 

In order to compare the model with our hom observations, we 
forward model the predicted emission from the loop simulation 
in each bandpass. This forward modeling is meant to provide 
a rough estimate of the relative change in emission created by 
embedding a single dynamic loop in a quiescent corona. First, 
we convolve the density and temperature distributions with the 
predicted AIA response function using the Solarsoft routine 
“aia_get_response.pro” (version 3), which specifically describes 
the spectral response of the instrument based on quiet-Sun level 
and ionization-equilibrium conditions and constant pressure. 
Next, we incorporate the resulting time-dependent intensity into 
a 3D structural model, which allows us to consider the projection 
effects of the surrounding ambient corona (Gibson et al. 2010). 
We assume the background emission is characterized by a 53° 
long cavity at a latitude of 60°, surrounded by a streamer with 
an isothermal temperature of 1.7 MK and a density profile 
taken from Schmit & Gibson (2011). The assumed extent of 
the evolving loop along the line of sight is 6 Mm. A study by 
Kucera et al. (2012) found that this approach overestimates the 
total photometric flux of coronal structures. 
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Figure 5. Time-distance diagram of TNE loop simulations. As in Figures 2 , (a) and (c) show 171 A emission in grayscale, while (b) and (d) show 193 A emission in 
grayscale. Regions strongly emitting in optically thin 304 A are shown in red hatching. Runs A and B differ in the magnitude of E\ (see Section 3.1). 

(A color version of this figure is available in the online journal.) 


Figure 5 shows the structures predicted by forward modeling 
the emission from simulation Runs A and B. We notice several 
similarities between the emission structures in these runs. While 
both runs exhibit brightenings in 171 A emission and similar 
condensation growth (Figures 5(a) and (c)), the intensity of 
the emission and the spatial extent varies; the change in 193 A 
emission (Figures 5(b) and (d)) is significantly weaker than the 
change in 17 1 A emission. Another important similarity between 
Runs A and B is in the temporal progression between the 
emission structures in different passbands. The 193 A emission 
brightens first, followed by an intensity increase in 171 A 
emission, followed by the formation of the condensation. This 
progression is consistent with earlier TNE studies (e.g., Luna 
et al. 2012). While the order of this progression is consistent, 
we find that the temporal offset between emissions in different 
bandpasses and the duration of emission strongly depends on 
the heating rate. 

Table 1 quantitatively compares model predictions of selected 
emission properties with their counterparts for observed promi- 
nence horns. One of the primary . differences between Run A and 
Run B is in the time lag between the peak emission in 171 A 
and the formation of the condensation, tm: stronger footpoint 
heating leads to a shorter time lag. For the eight-fold decrease 
in localized heating E\ between Runs A and B, we find that 
fi 7 i rises from 10 to 40 minutes. These brightenings exhibit dis- 
cernible rise and fall phases, which also depend on E \ : stronger 
heating yields shorter episodes of enhanced 171 A emission. The 
duration of the brightenings are related to the broad tempera- 
ture response functions for AIA and mass build up in the coronal 
portion of the loop. Run A accumulates mass in the corona more 
rapidly than Run B and therefore cools more rapidly through 
the peak of the temperature response function. Afni mea- 
sures the total duration of the brightening. The observationally 


Table 1 

Comparison of Model Values with Observations of Homs 



Run A 

Run B 

Homs 

log/ii 

-2.8 

-3.7 


/max (MK) 

3.4 

1.7 


*193 (hr) 

-0.4 

-2.9 

-0.5 ±0.4 

A *193 (hr) 

3.8 

2.2 

3 ±0.4 

A /193 

1.50 

1.01 

<1.2 

*171 (hr) 

0. 

-0.3 

-0.5 ±0.3 

A /171 (hr) 

2.1 

7.2 

3 ±0.3 

A /171 

11.0 

1.8 

1 . 2 - 1.8 


Notes. 7max is the maximum temperature reached during the simulation. *193 
is the time of peak emission in the 193 A bandpass, Afi 93 is the duration of the 
brightening, and A /193 is the maximum intensity of the brightening relative to 
background. All quantities subscripted “17 In” label congruent statistics, but for 
171 A emission. 

derived value of A*ni falls between the ranges provided by Runs 
A and B, while the observed value of *171 is longer than our 
simulations. For 193 A, the measured and simulated quantities 
agree well. 

One of the surprising attributes of the observed horns is the 
relative lack of correlation between the 171 A and 193 A emis- 
sion intensities, A /171 and A/ 193 . How could a thermodynamic 
perturbation create a 70% intensity perturbation in the 171 A 
bandpass and less than a 2% perturbation in the 193 A band- 
pass, given the large overlap in temperature response of these 
two bandpasses? The TNE model predicts that this discrepancy 
can be attributed to the nature of the cooling process. The tem- 
perature response functions for the 171 A and 193 A bandpasses 
each contain a single maximum within the coronal temperature 
range. The time of peak intensity for each bandpass occurs near, 
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but not at, that maximum emissivity temperature. Rather, the 
long-duration density increases in the corona shifts tn i and ^93 
to cooler temperatures with respect to their maximum emissiv- 
ity temperatures. The quantity A/193 is an order of magnitude 
smaller than Aim because the loop mass is steadily increasing 
while the local temperature drops from approximately 1.6 MK 
(peak emissivity in 193 A bandpass) to 0.8 MK (peak emissivity 
in 171 A bandpass). 

To summarize, we have compared the synthetic emission pre- 
dicted by the TNE model with the observed emission from 
prominence horns. Some model predictions match the obser- 
vations quite well, demonstrating that horns are a natural con- 
sequence of the condensation process treated in our model. In 
particular, the TNE model produces a 304 A emitting condensa- 
tion embedded within a coronal loop, and a 171 A brightening 
~50 Mm long and 3 hr in duration occurs 30 minutes before 
the formation of the condensation. The model also predicts that 
the 171 A brightening is much stronger than the 193 A bright- 
ening, as observed. Discrepancies also exist between the model 
and the observations, however, most notably in the progression 
of emission enhancement between bandpasses. We discuss the 
implications of this discrepancy in Section 5. 

4. CONSIDERATIONS FOR CAVITIES 

We now step back from the details of the comparisons 
described above and reassess the overarching question for 
the prominence-cavity system: what do the observations of 
prominence horns imply about the thermodynamic state of 
the cavity? In one school of thought, the cavity density is 
depleted because it feeds mass into the prominence (e.g., Saito 
& Tandberg-Hanssen 1973; Fort & Martres 1974; Berger et al. 
2012). This hypothesis is based on two distinct underlying 
assumptions: a condensation process occurs within the cavity 
and that process only converts mass from the coronal cavity 
into prominence plasma. The results of Section 3 suggest that 
horns are formed by cooling plasma, and horns appear to extend 
from the cavity into the prominence. Can this cooling process, 
when occurring over many loops, produce a large-scale density- 
depleted cavity? 

Tb directly answer this question, it is necessary to compare the 
plasma properties of the cavity with those of the predicted by our 
loop simulations. However, it is difficult to find common ground 
for a quantitative comparison due to key model assumptions. 
First, the loop geometry and the heating parameters applied in 
our model are suitable for studying the cooling and equilibration 
of a gravitationally stable condensation, whereas observed 
prominences and cavities typically are dynamic (Liu et al, 2012). 
Second, the condensation in the present simulations is produced 
by localized steady heating, whereas the actual temporal and 
spatial structure of coronal heating is a lively topic of debate 
(e.g., Klijnchuk 2006). The magnitude of the applied footpoint 
heating is within the range of canonical estimates (Withbroe & 
Noyes 1977), but this provides very loose constraints. In order 
to ultimately understand the time-averaged loop properties, 
we need more concrete information on the coronal heating 
mechanism itself and the lifecycle of prominence threads. 

Our simulations of TNE predict that condensation formation 
reduces the density in the coronal segment of the cooling 
loop, but only after the density has been increased through 
chromospheric evaporation. After a condensation has formed, 
the coronal portions of prominence-containing loops exhibit 
lower densities (70% at loop apex) and lower temperatures 
(50% at loop apex) than the Rosner et al. (1978) scaling laws 
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would predict for a half-loop of the same energy input. The 
TNE process can occur regardless of loop geometry, as long as 
the length-to-heating scale exceeds ~10 (Karpen et al. 2001). 
However, coronal rain is not prevalent in streamers, so the 
cavity-streamer density difference cannot be explained simply 
by invoking thermal non-equilibrium. One possibility is that 
different heating mechanisms operate in streamers and cavities, 
but thus far there is no compelling observational evidence for 
this suggestion. A change in magnetic structure also could 
differentiate between the cavity and the streamer. According 
to the Rosner et al. (1978) scaling laws, there is a weak direct 
relationship between loop length and density. These scaling 
laws do not apply to loops with the localized footpoint heating 
of the TNE model, however, or to dynamically evolving loops. 
Tb resolve this important issue, we need to determine whether 
the heating in the cavity differs significantly from that in 
the streamer, or whether geometrical effects can explain the 
difference. 

Better understanding of projection effects and filling factors 
also are important. Cavities can span hundreds of megameters 
in the axial direction. Given the dispersed distribution of 
high-latitude photospheric flux, a single large-scale cavity 
might encompass several distinct magnetic flux systems. In the 
synthetic emission shown in Section 3, the temporal- and spatial- 
averaged measurements were used to produce an emission 
background, but a more realistic approach would consider the 
filling factor of the condensing versus non-condensing field 
lines. Our prior multithread studies of prominences indicate that 
filling factors of order 0.001 yield reasonable agreement with 
estimated masses of observed prominences (Luna et al. 2012). 

One avenue of study that is worth pursuing is to probe for the 
presence of a relatively hot (>2 MK) component in the cavity. 
In most of the condensing loops studied here, the maximum 
loop temperature achieved prior to the cooling phase exceeds 
2 MK (Run A reaches 3.4 MK). This phase of the evolution 
may occur at a significantly lower density (depending on the 
equilibrium heating) than the cooling or post-condensation 
phase (see Figure 4), thus producing a weak EUV emission 
signature. The hot bandpasses available to AIA, specifically 
335 A and 131 A, exhibit low signal in the quiet Sun and the 
cavity. A long-exposure spectroscopic study using Hinode/ElS 
would allow us to set limits on the presence and variability of 
>2 MK plasma surrounding the prominence. 

This study suggests that horns are the manifestation of cooling 
plasma adding mass to the prominence. Based on the line-of- 
sight projection, these features appear inside the cavity. The 
TNE model predicts that a density depletion will occur in the 
corona that coincides with the formation of a condensation. 
While it is possible this depletion is responsible for the existence 
of the cavity, more evidence is needed to constrain the physical 
properties of the cavity and the streamer. 

5. DISCUSSION 

This analysis has combined two distinct but inherently related 
areas of prominence research: dynamic modeling and time- 
dependent EUV observations. This type of comparison is an 
essential step toward understanding the prominence-cavity 
system. We have attempted to present the predicted and observed 
data sets in a congruent manner and have found that our 
theoretical models are a close but inexact fit to the observations. 

We have presented evidence that, the dynamic emission 
observed surrounding prominences and projecting inside the 
cavity has several characteristics that are compatible with the 
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TNE model for a catastrophically cooling loop. The TNE 
provides relevant predictions on the thermodynamic state of the 
cavity, but a more thorough comparison requires a more detailed 
understanding of the heating mechanism and magnetic structure 
within and surrounding the prominence. The evidence at hand 
suggests that horns are related to the formation of prominence 
mass and that they connect the cavity and the prominence. 

One area of concern in our comparison of the TNE and the 
hom observations comes from the timing correlations between 
the various emission structures. Our simulations never exhibit 
171 A or 193 A emission that peaks after the formation of 
the condensation (such as Hom 3, Figure 2). One process 
that might produce this situation is a slowly cooling loop that 
undergoes a rapid mass increase coincidental with the formation 
of the condensation. Another possible scenario is that the line 
of sight intersects multiple loops in different stages of TNE, 
complicating the interpretation of the emission profiles from the 
horns, prominence, and cavity. Further study of our multithread 
prominence simulations could offer a more realistic approach 
to this 3D problem (e.g., Luna et al. 2012). 

Many avenues remain open to be explored in. our understand- 
ing of the prominence environment. The next step forward in 
our theoretical understanding is the self-consistent modeling of 
the magnetic structure and energetic evolution simultaneously, 
but this is a computationally challenging problem due to the 
large range of spatial and temporal scales that must be mod- 
eled. Moreover, while our ID loop simulations have expanded 
our understanding of the plasma dynamics and radiative signa- 
tures, the necessity of supplying such models with an empirical 
heating function will always limit the scope of their results. 
Choe & Lee (1992) offered a preliminary attempt at the full 
magnetic-plasma problem, namely that field line extension can 
siphon mass into the corona. Low et al. (2012) suggests that 
cross-field plasma transport may be important, highlighting an- 
other process that ID models cannot capture. It remains to be 
seen how the current generation of MHD models will elucidate 
the problem. The fundamental process described by TNE is 
that mass-loading the corona by evaporation of chromospheric 
plasma leads to cooling. Why, then, does this process occur 
in commonly in filament channels but infrequently in the sur- 
rounding streamer? Forthcoming observational studies of fila- 
ment channel and prominence evolution will shed light on the 
connection between magnetic structure and plasma dynamics, 
thus also illuminating the underlying physics of coronal heating 
and prominence formation. 
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